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PT-symmetric lattices with spatially extended gain/loss are 
generically unstable 
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PACS 42 . 25 . -p - Wave optics 

PACS 42 . 82 . Et - Waveguides, couplers, and arrays 

PACS 11.30.Er - Charge conjugation, parity, time reversal, and other discrete symmetries 

Abstract - We illustrate, through a series of prototypical examples, that linear parity-time (PT) 
symmetric lattices with extended gain/loss profiles are generically unstable, for any non-zero 
value of the gain/loss coefficient. Our examples include a parabolic real potential with a linear 
imaginary part and the cases of no real and constant or linear imaginary potentials. On the 
other hand, this instability can be avoided and the spectrum can be real for localized or compact 
PT-symmetric potentials. The linear lattices are analyzed through discrete Fourier transform 
techniques complemented by numerical computations. 
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Introduction. PT-symmetric quantum systems [T] (see 
also the review [2]) have emerged as an intriguing com- 
plex generalization of conventional quantum mechanics 
and have been a focus for numerous investigations in 
theoretical physics and applied mathematics. The key 
premise is that fundamental physical symmetries, such 
as parity P and time reversal T, may be sufficient (in 
suitable parametric regimes) to ensure that the eigenval- 
ues of the Hamiltonian are real. Thus, "PT-symmetric 
Hamiltonians provide an alternative to the standard pos- 
tulate that the Hamiltonian operator must be Hcrmi- 
tian. For Hamiltonians associated with a one-dimensional 
Schrodingcr operator with a complex potential V(x), the 
constraint of PT-symmctry requires that the potential 
satisfy V(x) = V(—x), i.e., V(x) has a symmetric real 
part and an antisymmetric imaginary part. 

Recently, the proposal that PT-symmetric systems can 
be physically implemented in the framework of optics [3] , 
was subsequently realized in experiments utilizing active 
or passive PT-dimers [3] and periodic lattices [5] . Similar 
proposals for the existence of a leaking dimcr (in the pres- 
ence of nonlinearity) have been formulated in the atomic 
setting of Bosc-Hubbard models [6]. Further experimen- 
tal investigations were concerned with electrical analogs 
of the linear PT-symmetric system [7] . Theoretical inves- 
tigations have rapidly followed by examining such dimer- 
type settings {5HTTJ] and generalizations thereof, including 



ones where the gain/loss contributions appear in a bal- 
anced form in front of the nonlinear term (11H13) . 

Although dimers have been the principal workhorse on 
which the investigation of an array of PT-symmetric sys- 
tems has been based, there are various works where more 
elements were considered [HI US]. PT-symmetric soli- 
tons were studied in full nonlinear lattices with a di- 
atomic structure [T6lJ18j . On the other hand, for the dimer 
monoatomic lattice, it is known that the PT-phase transi- 
tion occurs at gain/loss coefficient approaching zero when 
the number of lattice sites goes to infinity [19] . In other 
words, infinite lattices consisting of dimers are linearly un- 
stable and discrete solitons cannot be robust in such lat- 
tices, contrary to their counterparts in continuous models 
with periodic potentials [151120] . 

The purpose of the present work is to showcase the dra- 
matic differences between some prototypical discrete and 
continuum PT-symmetric systems. In particular, we show 
that linear PT-symmetric lattices with extended gain/loss 
are generically unstable: the eigenvalues of their associ- 
ated Hamiltonians are typically complex even if they are 
known to be purely real for Hamiltonians associated with 
the continuum analogue of the lattice. Our flagship ex- 
ample will be the standard quantum harmonic oscillator 
incorporating a purely imaginary linear potential [21j . Af- 
ter discretization, the spectrum of the discrete Schrodinger 
operator includes infinitely many complex eigenvalues, no 
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matter how small the gain/loss coefficient is. 

Other examples include the discrete Schrodinger oper- 
ator with no real potential, and constant or linear imag- 
inary potentials. For these examples, we also show that 
either infinitely many isolated eigenvalues or continuous 
spectral bands are complex inducing generic instability of 
linear 'PT-symmctric lattices in the presence of spatially 
extended gain/loss. 

Nevertheless, we also briefly touch upon a setting where 
an infinite lattice bearing PT-symmetry can have a real 
spectrum. This concerns the case of localized or compact 
potentials. For example, a PT-symmctric dimer embed- 
ded within the infinite lattice still enjoys real spectrum 
for sufficiently small gain/loss coefficient [22]. Therefore, 
localized or compact gain/loss may avoid the generic in- 
stability scenario, which we study here for lattices with 
spatially extended gain/loss. 

A potential with a parabolic real part and a linear imag- 
inary part. We start with the spectrum of the discrete 
Schrodinger operator with the potential V n = n 2 + i^n 
(cf. Rcf. [21] for the continuous version of this problem). 
Since V n = V- n , the potential is PT-symmetric. The rel- 
evant eigenvalue problem reads: 



Ellr, 



(u n+ i + u n _i - 2u n ) + (n 2 + i~/n) 



(1) 



Because the real part of the potential V n is bounded from 
below and confining asn->co, the spectrum of the linear 
lattice (p} is purely discrete [53] . The continuous version 
of the linear lattice Q]) takes the form: 



L = — 



dx 2 



„-7 
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Upon the change of variable x — > x + 17/2, the eigenvalue 
problem for L is tantamount to the quantum harmonic 
oscillator that has a purely real spectrum of eigenvalues 
located at E = 2m + 1 + 7 2 /4 for an integer m > [2"T] . 

On the contrary, as shown in Fig. [T] (top), the discrete 
case is significantly different in that the spectrum is pre- 
dominantly complex. The figure illustrates that, for differ- 
ent values of 7, there are only a few eigenvalues on the real 
axis and an infinite number of complex eigenvalues. The 
smaller the value of 7, the more eigenvalues are located 
on the real axis, but it is always a finite number for any 
7 7^ 0. The spectrum of the lattice ([1} is of course real for 
the Hermitian case of 7 = 0. 

To confirm these numerical findings with analytic the- 
ory, we introduce the discrete Fourier transform: 



1 r 

— / u{k)e- 



ikn 



dk. 



(2) 



Applying the discrete Fourier transform to the linear lat- 
tice ([!]) yields the differential equation in Fourier space: 



d 2 u dii r „ 
dk^ + ^dk + [E ' 



2cos(fc)]u(fc) = 0, 



(3) 
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Fig. 1: (Color online) Top panel: Eigenvalues of the linear lat- 
tice |T]), for the cases of 7 = 1 (blue circles), 7 = 0.5 (red stars) 
and 7 = 0.05 (green diamonds). The inset shows a blowup in 
the neighborhood of the smallest eigenvalues, highlighting the 
existence of one real eigenvalue for 7 = 1 and three real eigen- 
values for the other two cases. Bottom panel: trace of the 
monodromy matrix A(E) and the levels of 2cosh(-7r7) for the 
three values of 7. 



where we are looking for 27r-periodic functions u(k). Ap- 
plying the transformation v(k) = u(k)e lk / 2 , we obtain the 
Mathieu equation: 



cPv 
dk 2 



E - 2 - — + 2cos(fc) 
4 



v = 0. 



(4) 



Now we have v(k + 2ir) = e K1 v(k) : that is, we are looking 
for the Floquet multiplier /i» = e^ 7 of the monodromy 
matrix associated with the Mathieu equation (HJ). 

As it is well-known [24], the Floquet multiplier n(E) 
is determined from the trace of the monodromy matrix 
A(E) = fx(E) + fi(E)^ 1 . The function A diverges to pos- 
itive infinity as E — > — 00 and oscillates between values 
above 2 and below —2 for E > 0. Moreover, the local max- 
ima and minima of A(E) approach rapidly ±2 as E — > 00, 
because the cos(fc) potential in is smooth. 

The trace of monodromy matrix A(E) is shown in Fig.[T] 
(bottom) together with constant levels of 2 cosh(7T7) for 
the three values of 7. For any given 7 > 0, there are 
finitely many real roots E of equation /i(E) = /x* and the 
number of real roots grows as 7 — ► (it always includes 
at least one real root). For instance, we have one root for 
7=1 and three roots for 7 = 0.5 and 7 = 0.05. All other 
roots (infinitely many) are complex-valued. These roots 
bifurcate to complex numbers via saddle-node bifurcations 
when is increased from 1 (corresponding to 7 = 0) as 
7 is increased. 
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A potential with no real part and a linear imaginary 
part. Let us now drop the parabolic real potential and 
consider the linear eigenvalue problem: 
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Eu n = - (u n +l + lfn-1 ~ 2u n ) + i'jnUn- 



(5) 



Applying the discrete Fourier transform ([2]) yields the 
first-order differential equation in Fourier space: 



t4t + \E - 2 + 2 cos(fc)] u = 0. 
dk 

This equation can be exactly solved: 

u(k)=u(0)e^ 2 -V k - 2sin Wl 



(6) 



(7) 



The 27r-periodicity of the discrete Fourier transform u(k) 
gives now the eigenvalues E = 2 + i'ym, where m is an ar- 
bitrary integer. Hence, all the eigenvalues have an equidis- 
tant structure along the line Re(E) — 2 and again there 
are infinitely many complex eigenvalues. 

This conclusion can be checked directly from the differ- 
ence equation ([5]). If -Bo is an eigenvalue, then Eo + ijm is 
also an eigenvalue for any integer m thanks to the discrete 
group of symmetry of the linear lattice ([5]) with respect to 
translations in n. Also Eq = 2 is always an eigenvalue with 
the eigenvector u n = i n I n {2^~ 1 ), where /„ is the modified 
Besscl function. Note that the eigenvector is decaying in 
n because I n (2'y~ 1 ) — > as n —> ±oo for fixed 7 > 0. 

The sequence of equidistant eigenvalues along the line 
Re(-E) = 2 is illustrated in Fig. [2] for 7 = 1 (top) and 
7 = 0.1 (bottom). For 7 = 1, we find that the equidistant 
eigenvalues E = 2 + 17m are the only eigenvalues of the 
linear lattice ([5]). For smaller values of 7, we discover a 
new phenomenon, i.e., an appearance of additional parts 
of the spectrum of the linear lattice ([5]). For 7 = 0.1, the 
spectrum consists of the continuous spectrum that fills 
space between two curves to the right and left of the line 
Re(-E) = 2. For smaller values of 7, the relevant spec- 
tral "rectangle" keeps expanding along the real parts and 
contracting along the imaginary parts. The latter effect is 
induced by the finite-size truncation effects, because the 
spectrum of the infinite lattice always includes the equidis- 
tant eigenvalues E = 2 + ijm. The limit 7 — > is a sin- 
gular limit of the differential equation ©, when the only 
spectrum of the linear lattice (|5|) with 7 = is continuous 
and located on the real axis at [0,4]. 

We now prove that the sequence of equidistant eigenval- 
ues along the line Ke(E) = 2 is the only part of spectrum 
of the linear lattice (O for sufficiently large values of 7. 
Let E = 2 + 7A and rewrite the spectral problem ([5]) in 
the perturbed form: 



(A - in)u n = -7 1 (u n+ i + Un-i) 



(8) 



The limit 7—5-00 corresponds to the limit of weak cou- 
pling of the lattice and recovers the sequence of eigen- 
values at A = im for an integer m. Now, let A be any 
complex number different from {im} m ^z- Rewriting (j8J) 
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Fig. 2: (Color online) Eigenvalues of the linear lattice ([5}, for 
the cases of 7 = 1 (top) and 7 = 0.1 (bottom); the inset in the 
top panel shows the difference between the imaginary parts of 
adjacent eigenvalues) 

as u = r y~ 1 K(\)u, we can see that the operator K(\) 
is bounded and therefore, there exists a sufficiently large 
value of 7 such that this A cannot be in the spectrum of 
the linear lattice §5§. At the same time, A = im for any 
integer m is a simple eigenvalue for 7 = 00 that persists at 
A = im for any large 7. This analytical argument shows 
that no continuous spectrum exists for sufficiently large 7, 
so that there is a finite value of 7 = 70, for which a bi- 
furcation occurs in the spectrum of the linear lattice (JS|). 
Mathematical analysis of this bifurcation is beyond the 
scope of this work. 

A potential with no real part and a piecewise constant 
imaginary part. Let us now consider the piecewise con- 
stant potential in the linear eigenvalue problem: 



EUr, 



{u n+ i + u n -i - 2u„) + i7sign(n)u„ 



(9) 



To solve this linear problem, we introduce 9(E) as a 
complex- valued root of the dispersion relation 



E = 2 - 2cos((9) +17. 



(10) 



Since cos(#) is 2-7r-periodic and even, the root of the dis- 
persion relation (|10|) is uniquely determined in the semi- 
opened half-strip Re(#) € [— 7r,7r) and Im(#) > for any 
7 > and E e C with lm(E) ^ 7. Note that lm(0) = 
if lm(E) = 7. For any lm(E) 7^ 7, the eigenstate of the 
linear lattice ([9]) are represented by the exponentially de- 
caying function 



u n = 



e ie(E)n^ n > Q; 



(11) 



where 6*(E) is found from (|T0)) with the change 7 — > —7 
(and the same E). The linear lattice ((9]) is satisfied for 
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Fig. 3: (Color online) Eigenvalues of the linear lattice @, for 
the cases of 7 = 0.03 (top) and 7 = 0.06 (bottom). 

any n / 0, whereas the equation at n = gives another 
constraint on 9(E): 

e -ie(E)_ e -ie-(E) = iT (12) 

If Re(E) ^ 2, we find from (TO) and dHJ) that lm(6(E)) = 
lm(6(E)), which is impossible because they have opposite 
signs, hence Re(-E') = 2. Further studies of this system 
of equations with Re(E) = 2 show that the only solution 
exists for Im(E) = and corresponds to 

0(£O = ~+iarcsinQ|). (13) 

This construction gives the only eigenvalue of the linear 
lattice In addition, there are always two branches 

of the continuous spectrum located at lm(E) = ±7 for 
Re(E) £ [0,4], which correspond to the oscillatory non- 
decaying eigenfunctions with 6(E) £ R. We conclude that 
the only eigenvalue of the linear lattice (j9]) is stable, but 
the branches of the continuous spectrum are nevertheless 
unstable for any non-zero value of 7. 

Figure [3] shows eigenvalues of the truncated linear lat- 
tice © for 7 = 0.03 (top) and 7 = 0.06 (bottom). At a 
first glance, it seems that the numerical eigenvalues do not 
correspond to the analytical results above (one eigenvalue 
at E = 2 and the continuous spectrum at Im(E) = ±7 for 
Re(E) £ [0,4]). However, this is an artifact of the trunca- 
tion of the infinite lattice by a finite number TV of lattice 



sites subject to the Dirichlet end point conditions. When 
TV is increased from TV = 400 (stars) to TV = 800 (dots), 
the numerical eigenvalues are strongly affected and move 
towards the location of spectrum for the infinite lattice 
©■ 

Note that the spectrum of the linear lattice © is very 
different from the spectrum of the linear lattice of dimcrs: 

Eu n = - (lln+l + Un-i - 2u„) + ij(-l) n u n . (14) 

As it is well-known [T61ITH1I19] , the spectrum of the linear 
lattice of dimers (jT4")) is purely continuous and located at 



£ = 2±y4sin 2 i^-j - 7 2, 0e[-7T,7r], (15) 

which still shows instability for any non-zero value of 7 
because of the Fourier modes with 8 close to 0. 

A potential with a localized VT -symmetric part. Lastly, 
we mention briefly an example in which the PT-symmctric 
potential is localized in space, i.e., 

V n = scch 2 (n/M) (1 + i 7 tanh(n/M)) , (16) 

where parameter M denotes the width of the potential. 
If 7 = 0, the linear eigenvalue problem with the real po- 
tential (fT!)|) admits a continuous spectrum for E £ [0,4] 
and a finite number of simple eigenvalues for E > 4 
(because V n > 0, no eigenvalues appear for E < for 

7 = 0.) When 7 is small but non-zero, all simple eigenval- 
ues persist in PT-symmetric potentials by the perturba- 
tion theory similar to the one considered in |15j . On the 
other hand, since the potential V n decays exponentially 
as |n| — > 00, the continuous spectrum is not affected by 
the relatively compact perturbations and is determined 
by the oscillatory non-decaying eigenfunctions with real 

8 and E(6) =2 - 2cos(0) G [0,4]. Thus, we conclude 
that the spectrum of the linear eigenvalue problem with 
the PT-symmetric potential (Q~6|) is real at least for small 
7 > 0. 

Figure 2] shows the spectrum of the truncated linear 
lattice with the potential ([16]) for M = 10 and 7 = 0.02 
(top) or 7 = 0.1 (bottom). Simple eigenvalues for E > 4 
persist on the real axis for small values of 7. Although 
it seems that the continuous spectrum becomes complex 
near E = 0, this is again a numerical artifact, as adding 
more lattice sites TV reduces the imaginary part of the 
eigenvalues in the truncated lattice. As TV increases, we 
anticipate that the eigenvalues approach the real segment 
[0, 4] for the spectrum of the linear lattice with the poten- 
tial ([TBI). 

Note that the same conclusion holds also for the com- 
pact "PT-symmetric potentials such as the one correspond- 
ing to the embedded "PT-symmetric defect [22] : 

V n = ij(6 nt o - 5„ tl )u n , (17) 

where <5 n m is the Kronecker delta function. In this case, 
the eigenvalue spectrum spans again the interval [0,4] for 
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Fig. 4: (Color online) Eigenvalues of the linear lattice with 
the potential for M = 10 and 7 = 0.02 (top) or 7 = 0.1 
(bottom). 



all values of 7 e [—y/2, s/2\ in the infinite lattice size (i.e., 
N —> 00) limit. As an interesting aside, we mention that 
the critical point in this case (7 = y/2) is different from 
the critical point (of 7 = 1) for the PT-phase transition of 
a single dimer. In any case, we conclude that the compact 
support of the PT-symmetric defect enables a genuinely 
real spectrum at least for small values of gain/loss coeffi- 
cient in the infinite lattice limit, as was the case for the 
potential with the localized PT-symmetric part. 

Conclusions. We have considered the effect of discrete- 
ness on the spectral properties of linear PT-symmetric 
systems, characterized by spatially extended gain/loss and 
have contrasted these to the continuous PT-symmetric 
systems. We found that discreteness features generic in- 
stability of linear PT-symmetric lattices with extended 
gain/loss that can be unveiled by means of Fourier tech- 
niques, in conjunction with suitable analysis of the result- 
ing differential equations in Fourier space. We have also 
illustrated some remarkable differences in spectra between 
the finite lattice with N sites and the infinite lattice with 
N -> 00. 

The considered models with extended gain/loss featured 
generic instability properties that emerge essentially as 
soon as 7 7^ 0, i.e., their PT-phase transition occurs at 
7 = for all these linear lattices. Notice that the instabil- 
ity can not be avoided in the respective nonlinear settings, 



since PT-phase transition is a purely linear phenomenon, 
and nonlinearity can not arrest the exponential growth of 
the unstable eigenstates. Nevertheless, this instability can 
be dramatically avoided if the gain/loss are instead local- 
ized or compact: in this purely real spectrum of 
the infinite linear lattice can emerge for small values of 
gain/loss coefficient 7. 

There are many interesting directions for future stud- 
ies. First, providing a sharp criterion about the existence 
of real eigenvalues in infinite lattices depending on the lo- 
calization rate of the PT-symmetric potentials would be 
an extremely interesting condition both from a mathemat- 
ical and physical perspective. Furthermore, the difference 
in spectra between finite versus infinite lattices is another 
problem that merits additional investigation. It would also 
be interesting to extend these considerations to higher di- 
mensional lattices. 
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